------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /Users/sascha/Dropbox/suicide/RESTAT_replication/suicide-1816-restat-final.log
  log type:  text
 opened on:  22 Jul 2017, 10:11:26
r; t=0.00 10:11:26

. 
. set matsize 800
r; t=0.00 10:11:26

. set more off
r; t=0.00 10:11:26

. set rmsg on
r; t=0.00 10:11:26

. 
. 
. use "suicide-data-1816-restat.dta", clear
r; t=0.01 10:11:26

. 
. 
. global demog "f_pop1816_under15 f_pop1816_over60 f_urban1816"
r; t=0.00 10:11:26

. global econ  "buil1816_publ_pc"
r; t=0.00 10:11:26

. global educ  "edu1816_pri_enrol_i"
r; t=0.00 10:11:26

. global other "f_dead1816_1821_facc_all_pc"
r; t=0.00 10:11:26

. 
. /***************************************************/
. /* Table A.8: Descriptive Statistics, Prussia 1816 */
. /***************************************************/
. 
. # delimit ;
delimiter now ;
. global summvars "f_dead1816_1821_suic_all_pc f_dead1816_1821_suic_m_pc f_dead1816_1821_suic_f_pc f_dead1816_1821_suic_all_pd 
> f_rel1816_prot kmwittenberg $demog $econ $educ f_dead1816_1821_facc_all_pc f_dead1816_1821_facc_all_pd";
r; t=0.00 10:11:26

. # delimit cr
delimiter now cr
. 
. sum $summvars

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
f_~ic_all_pc |        306    6.501869    5.057208          0   26.05976
f_dead1~m_pc |        306    10.51911    8.280685          0   47.49842
f_dead1~f_pc |        306    2.690178    2.877741          0   22.42296
f_~ic_all_pd |        306     2.30211    1.899114          0   8.821367
f_rel1816_~t |        306    .5870899    .4096054          0   .9996294
-------------+---------------------------------------------------------
kmwittenberg |        306    .3160534    .1524554          0     .72962
f_pop1816~15 |        306     .361071     .028281   .2579834   .4557406
f_pop1816~60 |        306    .0658705    .0131714   .0389252   .1068461
 f_urban1816 |        306    .1198979    .2116877          0          1
buil1816_p~c |        306    .3278074    .3782693          0   2.090384
-------------+---------------------------------------------------------
edu1816_pr~i |        306    .5896957    .1974764   .0214847    1.10241
f_~cc_all_pc |        306    42.96322    16.72679   14.01095   123.5405
f_~cc_all_pd |        306     15.3918    6.765624   4.231312   54.89774
r; t=0.00 10:11:26

. 
. 
. /***********************************************************************/
. /* Table A.9: Protestantism and Suicide in Prussia 1816: OLS Estimates */
. /***********************************************************************/
. 
. regress f_dead1816_1821_suic_all_pc f_rel1816_prot, robust

Linear regression                               Number of obs     =        306
                                                F(1, 304)         =     183.31
                                                Prob > F          =     0.0000
                                                R-squared         =     0.3421
                                                Root MSE          =     4.1088

--------------------------------------------------------------------------------
               |               Robust
f_de~ic_all_pc |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
f_rel1816_prot |   7.221174   .5333543    13.54   0.000      6.17164    8.270707
         _cons |   2.262391   .2444099     9.26   0.000     1.781442     2.74334
--------------------------------------------------------------------------------
r; t=0.01 10:11:26

. regress f_dead1816_1821_suic_m_pc   f_rel1816_prot, robust

Linear regression                               Number of obs     =        306
                                                F(1, 304)         =     151.51
                                                Prob > F          =     0.0000
                                                R-squared         =     0.2996
                                                Root MSE          =     6.9416

--------------------------------------------------------------------------------
               |               Robust
f_dead181~m_pc |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
f_rel1816_prot |   11.06504    .898943    12.31   0.000     9.296104    12.83398
         _cons |   4.022937   .4155827     9.68   0.000     3.205154     4.84072
--------------------------------------------------------------------------------
r; t=0.00 10:11:26

. regress f_dead1816_1821_suic_f_pc   f_rel1816_prot, robust

Linear regression                               Number of obs     =        306
                                                F(1, 304)         =     123.61
                                                Prob > F          =     0.0000
                                                R-squared         =     0.2634
                                                Root MSE          =      2.474

--------------------------------------------------------------------------------
               |               Robust
f_dead181~f_pc |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
f_rel1816_prot |   3.605403   .3242859    11.12   0.000     2.967274    4.243533
         _cons |   .5734822   .1489974     3.85   0.000     .2802854     .866679
--------------------------------------------------------------------------------
r; t=0.00 10:11:26

. regress f_dead1816_1821_suic_all_pc f_rel1816_prot      $demog $econ $educ, robust

Linear regression                               Number of obs     =        306
                                                F(6, 299)         =      46.45
                                                Prob > F          =     0.0000
                                                R-squared         =     0.5094
                                                Root MSE          =     3.5777

-------------------------------------------------------------------------------------
                    |               Robust
f_dead181~ic_all_pc |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------------+----------------------------------------------------------------
     f_rel1816_prot |   4.734565   .5241358     9.03   0.000     3.703103    5.766028
  f_pop1816_under15 |  -36.22362   9.771961    -3.71   0.000    -55.45415   -16.99309
   f_pop1816_over60 |  -84.32181   18.09131    -4.66   0.000    -119.9242   -48.71939
        f_urban1816 |   5.977725    1.30222     4.59   0.000     3.415047    8.540402
   buil1816_publ_pc |   2.528509   .6590002     3.84   0.000     1.231642    3.825375
edu1816_pri_enrol_i |   3.135495   1.008774     3.11   0.002     1.150299    5.120691
              _cons |   18.96131   4.593029     4.13   0.000      9.92255    28.00007
-------------------------------------------------------------------------------------
r; t=0.00 10:11:26

. regress f_dead1816_1821_suic_m_pc   f_rel1816_prot      $demog $econ $educ, robust

Linear regression                               Number of obs     =        306
                                                F(6, 299)         =      38.58
                                                Prob > F          =     0.0000
                                                R-squared         =     0.4683
                                                Root MSE          =     6.0986

-------------------------------------------------------------------------------------
                    |               Robust
f_dead1816_182~m_pc |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------------+----------------------------------------------------------------
     f_rel1816_prot |   7.340539   .8785497     8.36   0.000     5.611615    9.069463
  f_pop1816_under15 |  -61.94906   16.80346    -3.69   0.000     -95.0171   -28.88102
   f_pop1816_over60 |  -134.6659   31.26408    -4.31   0.000    -196.1914   -73.14041
        f_urban1816 |   10.97393   2.451768     4.48   0.000     6.149019    15.79883
   buil1816_publ_pc |   2.487839   1.018922     2.44   0.015     .4826717    4.493007
edu1816_pri_enrol_i |   5.395588   1.759725     3.07   0.002     1.932573    8.858603
              _cons |   32.13504   7.980331     4.03   0.000     16.43031    47.83977
-------------------------------------------------------------------------------------
r; t=0.02 10:11:26

. regress f_dead1816_1821_suic_f_pc   f_rel1816_prot      $demog $econ $educ, robust

Linear regression                               Number of obs     =        306
                                                F(6, 299)         =      29.18
                                                Prob > F          =     0.0000
                                                R-squared         =     0.4049
                                                Root MSE          =     2.2422

-------------------------------------------------------------------------------------
                    |               Robust
f_dead1816_182~f_pc |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------------+----------------------------------------------------------------
     f_rel1816_prot |   2.232434   .3405458     6.56   0.000     1.562264    2.902605
  f_pop1816_under15 |  -12.44817    6.12126    -2.03   0.043    -24.49438   -.4019634
   f_pop1816_over60 |  -34.77439   12.07745    -2.88   0.004    -58.54196   -11.00682
        f_urban1816 |   1.679617    .814858     2.06   0.040     .0760337      3.2832
   buil1816_publ_pc |   2.584066    .761439     3.39   0.001     1.085608    4.082525
edu1816_pri_enrol_i |   1.069993   .6969911     1.54   0.126    -.3016365    2.441622
              _cons |   6.485392   2.893994     2.24   0.026     .7902152    12.18057
-------------------------------------------------------------------------------------
r; t=0.00 10:11:26

. 
. 
. /********************************************************************/
. /* Table 4: Protestantism and Suicide in Prussia 1816: IV Estimates */
. /********************************************************************/
. 
. /* IV 1st stage */
. regress f_rel1816_prot  kmwittenberg  $demog $econ $educ, robust

Linear regression                               Number of obs     =        306
                                                F(6, 299)         =      64.61
                                                Prob > F          =     0.0000
                                                R-squared         =     0.3603
                                                Root MSE          =     .33089

-------------------------------------------------------------------------------------
                    |               Robust
     f_rel1816_prot |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------------+----------------------------------------------------------------
       kmwittenberg |  -.8688479    .195125    -4.45   0.000     -1.25284   -.4848557
  f_pop1816_under15 |  -4.250666    .942651    -4.51   0.000    -6.105737   -2.395595
   f_pop1816_over60 |  -8.089395   1.844953    -4.38   0.000    -11.72013   -4.458656
        f_urban1816 |  -.0781313   .0876583    -0.89   0.373    -.2506366     .094374
   buil1816_publ_pc |   .0914158   .0563807     1.62   0.106    -.0195375    .2023691
edu1816_pri_enrol_i |   .4433432   .1040283     4.26   0.000     .2386229    .6480636
              _cons |   2.647301   .4207934     6.29   0.000     1.819209    3.475392
-------------------------------------------------------------------------------------
r; t=0.01 10:11:26

. 
. /* IV 2nd stage */
. ivreg2 f_dead1816_1821_suic_all_pc  (f_rel1816_prot=kmwittenberg) $demog $econ $educ, robust

IV (2SLS) estimation
--------------------

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity

                                                      Number of obs =      306
                                                      F(  6,   299) =    28.58
                                                      Prob > F      =   0.0000
Total (centered) SS     =  7800.483711                Centered R2   =   0.0244
Total (uncentered) SS   =  20736.41924                Uncentered R2 =   0.6330
Residual SS             =  7610.183141                Root MSE      =    4.987

-------------------------------------------------------------------------------------
                    |               Robust
f_dead181~ic_all_pc |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------------+----------------------------------------------------------------
     f_rel1816_prot |   14.98912   3.020265     4.96   0.000     9.069506    20.90873
  f_pop1816_under15 |   15.25326   21.53809     0.71   0.479    -26.96062    57.46715
   f_pop1816_over60 |   1.280613   39.49954     0.03   0.974    -76.13706    78.69828
        f_urban1816 |   6.570462    1.43104     4.59   0.000     3.765676    9.375248
   buil1816_publ_pc |  -.2621514   1.165411    -0.22   0.822    -2.546315    2.022012
edu1816_pri_enrol_i |   -2.07355   1.997866    -1.04   0.299    -5.989294    1.842195
              _cons |  -7.369042   10.83011    -0.68   0.496    -28.59567    13.85758
-------------------------------------------------------------------------------------
Underidentification test (Kleibergen-Paap rk LM statistic):             26.404
                                                   Chi-sq(1) P-val =    0.0000
------------------------------------------------------------------------------
Weak identification test (Cragg-Donald Wald F statistic):               29.577
                         (Kleibergen-Paap rk Wald F statistic):         19.827
Stock-Yogo weak ID test critical values: 10% maximal IV size             16.38
                                         15% maximal IV size              8.96
                                         20% maximal IV size              6.66
                                         25% maximal IV size              5.53
Source: Stock-Yogo (2005).  Reproduced by permission.
NB: Critical values are for Cragg-Donald F statistic and i.i.d. errors.
------------------------------------------------------------------------------
Hansen J statistic (overidentification test of all instruments):         0.000
                                                 (equation exactly identified)
------------------------------------------------------------------------------
Instrumented:         f_rel1816_prot
Included instruments: f_pop1816_under15 f_pop1816_over60 f_urban1816
                      buil1816_publ_pc edu1816_pri_enrol_i
Excluded instruments: kmwittenberg
------------------------------------------------------------------------------
r; t=0.06 10:11:26

. ivreg2 f_dead1816_1821_suic_all_pc  (f_rel1816_prot=kmwittenberg) $demog $econ $educ $other, robust

IV (2SLS) estimation
--------------------

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity

                                                      Number of obs =      306
                                                      F(  7,   298) =    24.92
                                                      Prob > F      =   0.0000
Total (centered) SS     =  7800.483711                Centered R2   =   0.0261
Total (uncentered) SS   =  20736.41924                Uncentered R2 =   0.6337
Residual SS             =  7596.502384                Root MSE      =    4.982

---------------------------------------------------------------------------------------------
                            |               Robust
f_dead1816_1821_suic_all_pc |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
----------------------------+----------------------------------------------------------------
             f_rel1816_prot |   14.97072   2.873271     5.21   0.000     9.339216    20.60223
          f_pop1816_under15 |   15.11707   20.66683     0.73   0.464    -25.38917     55.6233
           f_pop1816_over60 |   .5549999   36.46381     0.02   0.988    -70.91276    72.02276
                f_urban1816 |   6.565985   1.427905     4.60   0.000     3.767342    9.364628
           buil1816_publ_pc |  -.2562612   1.128986    -0.23   0.820    -2.469033    1.956511
        edu1816_pri_enrol_i |  -2.044539    1.88103    -1.09   0.277    -5.731291    1.642213
f_dead1816_1821_facc_all_pc |  -.0012826   .0205212    -0.06   0.950    -.0415035    .0389383
                      _cons |  -7.224666   10.02869    -0.72   0.471    -26.88053     12.4312
---------------------------------------------------------------------------------------------
Underidentification test (Kleibergen-Paap rk LM statistic):             31.418
                                                   Chi-sq(1) P-val =    0.0000
------------------------------------------------------------------------------
Weak identification test (Cragg-Donald Wald F statistic):               33.991
                         (Kleibergen-Paap rk Wald F statistic):         24.861
Stock-Yogo weak ID test critical values: 10% maximal IV size             16.38
                                         15% maximal IV size              8.96
                                         20% maximal IV size              6.66
                                         25% maximal IV size              5.53
Source: Stock-Yogo (2005).  Reproduced by permission.
NB: Critical values are for Cragg-Donald F statistic and i.i.d. errors.
------------------------------------------------------------------------------
Hansen J statistic (overidentification test of all instruments):         0.000
                                                 (equation exactly identified)
------------------------------------------------------------------------------
Instrumented:         f_rel1816_prot
Included instruments: f_pop1816_under15 f_pop1816_over60 f_urban1816
                      buil1816_publ_pc edu1816_pri_enrol_i
                      f_dead1816_1821_facc_all_pc
Excluded instruments: kmwittenberg
------------------------------------------------------------------------------
r; t=0.05 10:11:26

. ivreg2 f_dead1816_1821_suic_m_pc    (f_rel1816_prot=kmwittenberg) $demog $econ $educ $other, robust

IV (2SLS) estimation
--------------------

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity

                                                      Number of obs =      306
                                                      F(  7,   298) =    22.92
                                                      Prob > F      =   0.0000
Total (centered) SS     =   20913.7727                Centered R2   =   0.0228
Total (uncentered) SS   =  54773.19047                Uncentered R2 =   0.6269
Residual SS             =   20436.0128                Root MSE      =    8.172

---------------------------------------------------------------------------------------------
                            |               Robust
  f_dead1816_1821_suic_m_pc |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
----------------------------+----------------------------------------------------------------
             f_rel1816_prot |   23.43898   4.701312     4.99   0.000     14.22458    32.65339
          f_pop1816_under15 |   19.21804   33.92617     0.57   0.571    -47.27603    85.71211
           f_pop1816_over60 |   4.341455   59.44552     0.07   0.942    -112.1696    120.8525
                f_urban1816 |   11.93203   2.536181     4.70   0.000     6.961207    16.90285
           buil1816_publ_pc |   -1.90032   1.706037    -1.11   0.265    -5.244092    1.443452
        edu1816_pri_enrol_i |   -2.94089   3.160456    -0.93   0.352    -9.135271     3.25349
f_dead1816_1821_facc_all_pc |   .0103617     .03361     0.31   0.758    -.0555126    .0762361
                      _cons |   -9.98536   16.50173    -0.61   0.545    -42.32816    22.35744
---------------------------------------------------------------------------------------------
Underidentification test (Kleibergen-Paap rk LM statistic):             31.418
                                                   Chi-sq(1) P-val =    0.0000
------------------------------------------------------------------------------
Weak identification test (Cragg-Donald Wald F statistic):               33.991
                         (Kleibergen-Paap rk Wald F statistic):         24.861
Stock-Yogo weak ID test critical values: 10% maximal IV size             16.38
                                         15% maximal IV size              8.96
                                         20% maximal IV size              6.66
                                         25% maximal IV size              5.53
Source: Stock-Yogo (2005).  Reproduced by permission.
NB: Critical values are for Cragg-Donald F statistic and i.i.d. errors.
------------------------------------------------------------------------------
Hansen J statistic (overidentification test of all instruments):         0.000
                                                 (equation exactly identified)
------------------------------------------------------------------------------
Instrumented:         f_rel1816_prot
Included instruments: f_pop1816_under15 f_pop1816_over60 f_urban1816
                      buil1816_publ_pc edu1816_pri_enrol_i
                      f_dead1816_1821_facc_all_pc
Excluded instruments: kmwittenberg
------------------------------------------------------------------------------
r; t=0.05 10:11:26

. ivreg2 f_dead1816_1821_suic_f_pc    (f_rel1816_prot=kmwittenberg) $demog $econ $educ $other, robust

IV (2SLS) estimation
--------------------

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity

                                                      Number of obs =      306
                                                      F(  7,   298) =    16.77
                                                      Prob > F      =   0.0000
Total (centered) SS     =  2525.824584                Centered R2   =   0.0777
Total (uncentered) SS   =  4740.364334                Uncentered R2 =   0.5086
Residual SS             =  2329.559817                Root MSE      =    2.759

---------------------------------------------------------------------------------------------
                            |               Robust
  f_dead1816_1821_suic_f_pc |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
----------------------------+----------------------------------------------------------------
             f_rel1816_prot |   7.066208    1.62077     4.36   0.000     3.889556    10.24286
          f_pop1816_under15 |   11.34511   11.62027     0.98   0.329    -11.43021    34.12042
           f_pop1816_over60 |  -.5757935   21.74905    -0.03   0.979    -43.20314    42.05155
                f_urban1816 |   1.922307   .9001465     2.14   0.033     .1580526    3.686562
           buil1816_publ_pc |   1.278127   1.032823     1.24   0.216    -.7461697    3.302424
        edu1816_pri_enrol_i |  -1.173917   .9631684    -1.22   0.223    -3.061692    .7138585
f_dead1816_1821_facc_all_pc |  -.0137942   .0111547    -1.24   0.216    -.0356571    .0080687
                      _cons |  -4.881345   5.699649    -0.86   0.392    -16.05245    6.289761
---------------------------------------------------------------------------------------------
Underidentification test (Kleibergen-Paap rk LM statistic):             31.418
                                                   Chi-sq(1) P-val =    0.0000
------------------------------------------------------------------------------
Weak identification test (Cragg-Donald Wald F statistic):               33.991
                         (Kleibergen-Paap rk Wald F statistic):         24.861
Stock-Yogo weak ID test critical values: 10% maximal IV size             16.38
                                         15% maximal IV size              8.96
                                         20% maximal IV size              6.66
                                         25% maximal IV size              5.53
Source: Stock-Yogo (2005).  Reproduced by permission.
NB: Critical values are for Cragg-Donald F statistic and i.i.d. errors.
------------------------------------------------------------------------------
Hansen J statistic (overidentification test of all instruments):         0.000
                                                 (equation exactly identified)
------------------------------------------------------------------------------
Instrumented:         f_rel1816_prot
Included instruments: f_pop1816_under15 f_pop1816_over60 f_urban1816
                      buil1816_publ_pc edu1816_pri_enrol_i
                      f_dead1816_1821_facc_all_pc
Excluded instruments: kmwittenberg
------------------------------------------------------------------------------
r; t=0.18 10:11:26

. ivreg2 f_dead1816_1821_suic_all_pd  (f_rel1816_prot=kmwittenberg) $demog $econ $educ $other, robust

IV (2SLS) estimation
--------------------

Estimates efficient for homoskedasticity only
Statistics robust to heteroskedasticity

                                                      Number of obs =      306
                                                      F(  7,   298) =    26.10
                                                      Prob > F      =   0.0000
Total (centered) SS     =  1100.023812                Centered R2   =   0.0921
Total (uncentered) SS   =  2721.735044                Uncentered R2 =   0.6331
Residual SS             =  998.6709822                Root MSE      =    1.807

---------------------------------------------------------------------------------------------
                            |               Robust
f_dead1816_1821_suic_all_pd |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
----------------------------+----------------------------------------------------------------
             f_rel1816_prot |   5.239947   1.059419     4.95   0.000     3.163523    7.316371
          f_pop1816_under15 |   4.804406   7.797329     0.62   0.538    -10.47808    20.08689
           f_pop1816_over60 |   7.929221    13.5925     0.58   0.560     -18.7116    34.57004
                f_urban1816 |   1.809194   .4932976     3.67   0.000      .842348    2.776039
           buil1816_publ_pc |   .0872183   .3944288     0.22   0.825    -.6858478    .8602845
        edu1816_pri_enrol_i |  -.8861152   .6762981    -1.31   0.190    -2.211635    .4394048
f_dead1816_1821_facc_all_pc |   .0001464   .0083721     0.02   0.986    -.0162627    .0165555
                      _cons |  -2.760506   3.755163    -0.74   0.462    -10.12049    4.599478
---------------------------------------------------------------------------------------------
Underidentification test (Kleibergen-Paap rk LM statistic):             31.418
                                                   Chi-sq(1) P-val =    0.0000
------------------------------------------------------------------------------
Weak identification test (Cragg-Donald Wald F statistic):               33.991
                         (Kleibergen-Paap rk Wald F statistic):         24.861
Stock-Yogo weak ID test critical values: 10% maximal IV size             16.38
                                         15% maximal IV size              8.96
                                         20% maximal IV size              6.66
                                         25% maximal IV size              5.53
Source: Stock-Yogo (2005).  Reproduced by permission.
NB: Critical values are for Cragg-Donald F statistic and i.i.d. errors.
------------------------------------------------------------------------------
Hansen J statistic (overidentification test of all instruments):         0.000
                                                 (equation exactly identified)
------------------------------------------------------------------------------
Instrumented:         f_rel1816_prot
Included instruments: f_pop1816_under15 f_pop1816_over60 f_urban1816
                      buil1816_publ_pc edu1816_pri_enrol_i
                      f_dead1816_1821_facc_all_pc
Excluded instruments: kmwittenberg
------------------------------------------------------------------------------
r; t=0.05 10:11:26

. 
. 
. /**********************************************************************************************************/
. /**********************************************************************************************************/
. 
. 
. 
. log close
      name:  <unnamed>
       log:  /Users/sascha/Dropbox/suicide/RESTAT_replication/suicide-1816-restat-final.log
  log type:  text
 closed on:  22 Jul 2017, 10:11:26
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
